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Abstract 

Motivated by a phenomenon of phase transition in a model of ahgnment 
of self-propehed particles, we obtain a kinetic mean-field equation which is 
nothing else than the Doi equation (also called Smoluchowski equation) with 
dipolar potential. 

In a self-contained article, using only basic tools, we analyze the dynamics 
of this equation in any dimension. We first prove global well-posedness of this 
equation, starting with an initial condition in any Sobolev space. We then 
compute all possible steady-states. There is a threshold for the noise parame- 
ter: over this threshold, the only equilibrium is the uniform distribution, and 
under this threshold, there is also a family of non-isotropic equilibria. 

We give a rigorous prove of convergence of the solution to a steady-state 
as time goes to infinity. In particular we show that in the supercritical case, 
the only initial conditions leading to the uniform distribution in large time are 
those with vanishing momentum. For any positive value of the noise parame- 
ter, and any initial condition, we give rates of convergence towards equilibrium, 
exponentially for both supercritical and subcritical cases and algebraically for 
the critical case. 
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1 Introduction 



Phase transition and large time behavior of large interacting oriented/rod-like parti- 
cle systems and their mean field limits have shown to be interesting in many physical 
and biological complex systems. Examples are: paramagnetism to ferromagnetism 
phase transition near Curie temperature, nematic phase transition in liquid crys- 
tal or rod-shaped polymers, emerging of fiocking dynamics near critical mass of 
self-propelled particles, etc. 

The dynamics on orientation for self-propelled particles proposed by Vicsek et 
al [12] to describe, for instance, fish schooling or bird fiocking, present such a behav- 
ior in numerical simulations. As the density increases (or as the noise decreases) and 
reaches a threshold one can observe strong correlations between the orientations of 
particles. The model is discrete in time and particles move at constant speed follow- 
ing their orientation. At each time step, the orientation of each particle is updated, 
replaced by the mean orientation of its neighbors, plus a noise term. 

A way to provide a time-continuous version of this dynamical system, which 
allows to take a mean-field limit (and even a macroscopic limit), has been proposed 
by Degond and Motsch [6]. Instead of replacing the orientation at the next time 
step, they introduce a parameter playing the role of a rate of relaxation towards this 
mean orientation. Unfortunately the mean-field limit of this model does not present 
phase transition. In [TT], the first author of the present paper proved the robustness 
of the behavior of this model when this rate of relaxation depends on a local density. 
In particular, phase transition is still absent. However, when this parameter is set 
to be proportional to the local momentum of the neighboring particles, we will see 
that the model present a phenomenon of phase transition. This phenomenon occurs 
on the orientation dynamics, so we will only consider here the spatial homogeneous 
dynamics. A current joint work with Pierre Degond is dedicated to the study of this 
model when we take in account the space variable. It is left aside in the present 
paper, and a specific paper j5] is in progress on this subject. 

The particular model is described as follows: we have oriented particles, de- 
scribed by vectors Wi, . . . jOJat belonging to S, the unit sphere of M", and satisfying 
the following system of coupled stochastic differential equations, for k G |1, A^J: 



The term (Id — uok® ujk) denotes the projection on the hyperplane orthogonal to Uk, 
and constrains the norm of Uk to be constant. The terms stand for A^ independent 
standard Brownian motions on M", and then the stochastic term (Id — Uk ®Uk)dB^ 
represents the contribution of a Brownian motion on the sphere S to the model. For 
more details on how to define Brownian motion on a Riemannian manifold, see |12] . 
Without this stochastic term, equation ([T]) can be written 



where V^; is the tangential gradient on the sphere (see the beginning of Section 12.11 
for some useful formulas on the unit sphere). So the model can be understood as a 






(2) 
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relaxation towards a unit vector in the direction of J^, subjected to a Brownian mo- 
tion on the sphere with intensity -\/2a. The only difference with the model proposed 
in [6] (in the spatial homogeneous case) is that Jk is there replaced by vVLk-, where Vtk 
is the unit vector in the direction of and the frequency of relaxation p is constant 
(or dependent on the local density in [H]). One point to emphasize is that, in that 
model, the interaction cannot be seen as a sum of binary interactions, contrary to 
the model presented here. Here the mean momentum Jk does not depend on the 
index k (but this is not true in the inhomogeneous case, where the mean is taken 
among the neighboring particles). 

To simplify notations, we work with the uniform measure of total mass 1 on the 
sphere §. We denote by : M+ x S — )■ the probability density function (de- 
pending on time) associated to the position of one particle. Then, as the number 
of particles tends to infinity, tends to a probability density function / satisfying 

dtf = QU). (3) 

with 

QU) = -V. ■ ((Id -u® uj)J[f]f) + aA^/, (4) 

J[f]= ! ujf{.,uj)duj. (5) 

In the model of |6], J[f] is just replaced in (jl]) by z/f2[/], where ^[f] is the unit 
vector in the direction of J[f]- 

The first term of Q{f) can be formally derived using a direct computation with 
the empirical distribution of particles. And the diffusion part comes from Ito's 
formula. A rigorous derivation of this mean-field limit is outside the scope of the 
present paper, and is linked with the so-called "propagation of chaos" property. We 
refer to [18] for an introduction to this notion. The laboratory example given in this 
reference is the original model of McKean [T3] which is a more general version of 
our system in instead of § (in that case, equation ([H]) is called McKean- Vlasov 
equation). The main point is to adapt the theory in the framework of stochastic 
analysis on Riemannian manifolds. 

Notice that equation ([3]) can be written in the form 

dtf = V- (/Vvl>) + aA/, 

with 

^{u,t) = -u ■ J{t) = K{uj,u}) f{t,u})du!. 

Js 

This equation is known as Doi equation (or Doi-Onsager, Smoluchowski, or even 
nonlinear Fokker-Planck equation) and was introduced by Doi [7] as gradient flow 
equation for the Onsager free energy functional: 

nf) = a [ f{.,u)\nf{.,u)du + l [ K{u,u)f{.,u)f{.,u)dudu. (6) 

This functional was proposed by Onsager to describe the equilibrium states 
of suspensions of rod-like polymers. They are given by the critical points of this 
functional. 
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Defining the chemical potential fi as the first order variation of J^{f) under the 
constraint Jg / = 1, we get /i = cxln / + and the Doi equation becomes 



In the original work of Onsager, the kernel has the form K{uj,u}) = \u x uj\, 
but there is another form, introduced later by Maier and Saupe [13], which leads to 
similar quantitative results: K{uj,oj) = —{u ■ a))^. In our case, the potential given 
by K{oj,u}) = —oj ■ Co is called the dipolar potential. This is a case where the arrow 
of the orientational direction has to be taken in account. 

One of the interesting behavior of the Doi-Onsager equation is the phase tran- 
sition bifurcation. This is indeed easy to see (here with the dipolar potential) from 
the following linearization around the uniform distribution: if / is a probability 
density function, solution of (l3l), we write f = 1 + g, so j^gdu = and we can 
get the equation for g. We multiply the equation by uj and integrate, using the 
formula J^u ^ u do; = i Id (this is a matrix with trace one and commuting with 
any rotation) and the tools in the beginning of Section 12.11 We get the linearized 
equation for g and J[g]: 



Therefore if we take the linear part of this system, we can solve the second 
equation directly, and the first one becomes the heat equation with a known source 
term. Finally, around the constant state, the linearized Doi equation is stable if 
a ^ ^, and unstable if cr < j^. We expect to find another kind of equilibrium in this 
regime. The work has been done in jH] for the dimension n = 3, the distribution 
obtained is called Fisher- Von Mises distribution [20] • As far as we know, this is the 
only work dealing with the dipolar potential alone. 

A lot of work has been done to study the equilibrium states for the Maier-Saupe 
potential, and in particular to show the axial symmetry of these steady states. A 
complete classification has been achieved for the two and three-dimensional cases 
in [13] (see also [21], including the analysis of stability under a weak external shear 
flow). The interesting behavior, besides the phase transition, is the hysteresis phe- 
nomenon: before a first threshold, only the anisotropic equilibrium is stable, then 
both anisotropic and uniform equilibria are stable, and after a second threshold, 
the only equilibrium is the uniform distribution. In the case of a coupling between 
the Maier-Saupe and the dipolar, it is shown in [23] that the only stable equilib- 
rium states are axially symmetric. To our knowledge, less work has been done to 
study the dynamics of the Doi-Onsager equation, in particular the rate at which the 
solution converges to a steady-state. 

The purpose of this paper is to give a rigorous proof of the phase transition in 
any dimension for the dipolar potential, and study the large time dynamics and the 
convergence rates towards equilibrium states. 

In Section [21 we give some general results concerning equation ([3]). We provide 
a self-contained proof for existence and uniqueness of a solution with initial nonneg- 
ative condition in any Sobolev space. We show that the solution is instantaneously 



5*/ = V-(/V/i). 



dtg = (xA^g +{n-l)u- J[g] + 0{g^), 
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positive and in any Sobolev space (and actually analytic in the space variable), and 
we obtain uniform bounds in time for each Sobolev norm. 

In Section [3l we use the Onsager free energy (decreasing in time) to analyze the 
general behavior of the solution as time goes to infinity. We prove a kind of LaSalle 
principle, implying that the solution converges, in the cj-limit sense, to a given set 
of equilibria. We determine all the steady states, and see that the value ^ is indeed 
a threshold for the noise parameter a. Over this threshold, the only equilibrium is 
the uniform distribution. When cr < two kinds of equilibria exist: the uniform 
distribution, and a family of non-isotropic distributions (called Fischer- Von Mises 
distributions), with a concentration parameter k depending on a. 

Finally, in Section |H we show that the solution converges strongly to a given 
equilibrium. We first obtain a new conservation relation, which plays the role of an 
entropy when cr ^ -, and shows a global convergence to the uniform distribution 
with rate proportional to o—^. Then we prove that, in the supercritical case cr < ^, 
the solution converges to a non-isotropic equilibrium if and only if the initial drift 
velocity |</[/o]| is non-zero (if it is zero, the equation reduces to the heat equation, 
and the solution converges exponentially fast to the uniform distribution). We prove 
in that case that the convergence to this steady-state is exponential in time, and 
we give the asymptotic rate of convergence. Finally, in the critical case cr = ^, we 
show that the speed of convergence to the uniform distribution is algebraic (more 
precisely the decay in any Sobolev norm is at least ;^). 



2 General results 



2.1 Preliminaries: some results on the unit sphere 

This subsection consists essentially in a main lemma, allowing to perform some 
estimates on the norm of integrals of the form J^gV^jh, where h and g are real 
functions with mean zero. 

But let us start by some useful formulas. 

For V a constant vector in M", we have: 

V^(a;- V) = (Id-uj^uj)V 
■ ((Id - a; ® u)V) = -{n-l)uj-V, 

where V^^ (resp. V^j-) stands for the tangential gradient (resp. the divergence) on 
the unit sphere. When no confusion is possible, we will just use the notation V. 

Then, taking the dot product with a given tangent vector field A or multiplying 
by a regular function / and integrating by parts, we get 

ujVuj- A{uj)du = - A{uj)du 

I VJduj = {n-l) [ ufdu. 

We then introduce some notations. We denote by if''(§) the subspace composed 
of mean zero functions of the Sobolev space i^'^(§). This is a Hilbert space, associ- 
ated to the inner product {g, ^) = {i~'^)^9> where A is the Laplace-Beltrami 
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operator on the sphere. This has also a sense for any s G M by spectral decomposition 
of this operator. We will denote by || ■ the norm on this Hilbert space. 

We then define the so-called conformal Laplacian A„_i on the sphere (see jl]) 
which plays a role in some Sobolev inequalities. This is a positive definite op- 
erator (pseudodifferential operator of degree n — 1, mapping continuously H'^(S) 
into which is a differential operator when n is odd) given by 



H {-A+j{n-j-2)) 



for n odd. 



(-A + (f - 1)^) ^ Yl (-A + j{n -j- 2)) for n even. 



(7) 



Equivalently, it can be also defined by 

An_i = £{i + 1) ...(£ + n — 2)1^ for any spherical harmonic of degree i. (8) 

Here is the main lemma. 

Lemma 1. Estimates on the sphere. 

1. If h in H^^^^(§>) and g in H'^(S>), the following integral is well defined and we 
have 



gVh 



^ C'll£'lli/-ll^lli/-+i 



(9) 



where the constant C depends only on s and n. 
2. We have the following estimation, for any g G H^^^{^): 



[ gVi-AYg 



^C\\g\\l^ 



(10) 



where the constant C depends only on s and n. 
3. We have the following identity, for any g G 

gVA-\g = 



Let us make some remarks on these statements. The first one is just expressing 
the fact that the gradient operator (or more precisely any of its component e • V for 
a given unit vector e) is well defined as an operator sending i/^*+^(S) continuously 
into H~^{E>) for any s. 

The second one is actually a commutator estimate. It is equivalent to the fact 
that for any given unit vector e, and for any g,h E H^^^ we have 



ge-V{-Ayh + he-V{-Ayg 



^c\\gum 



Defining the operator F by 

Fg = e- V{-Ayg - (-A)^V • ((Id - w ® io)eg) 
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and integrating by parts, this inequality becomes \J^hFg\ ^ In 
other words, F sends H''^{S) continuously into H~^{S) for any s. 

So since F = [e ■ V, (—A)''] + (n — 1)(— A)'^e ■ u, this second statement (fTOl) 
expresses that the commutator [V, (—A)''] is an operator of degree 2s. 

With the same point of view, the last equality fill I) gives an exact computation 
of the commutator of the gradient and the inverse of conformal Laplacian. 

This is just saying that [V, A"^;^] = — (n — l)A~\u}, or, multiplying left and 
right by A„_i, that [V, A„„i] = (n - l)a;A„_i. 

The proof of this lemma relies on some computations on spherical harmonics, 
and is given in Appendix lA.ll 

2.2 Existence, uniqueness, positivity, regularity. 

We present here a self-contained proof of well-posedness of the problem ([3]), working 
in any Sobolev space for the initial condition. Some analogous claims are given in 
without proof, starting for a continuous nonnegative function. They are based on 
arguments of [2], stating that the Galerkin method based on spherical harmonics 
converges (exponentially fast) to the unique solution. They are weaker with respect 
to the initial conditions and the positivity, but stronger for the regularity of the 
solution (analytic in space). As a remark we will give the same regularity results, 
and prove it in Appendix IA.2[ 

Definition 2.1. Weak solution. 

ForT > 0, the function f G L\{0,T), H'+\S)) H H\{0,T), H'-\S)) is said to 
be a weak solution of ([3]) if for almost all t G [0, T], we have for all h G H~^~^^(E>) 

{dj,h) = -a{V^f,V^h) + (/, J[/] ■ V^/i), (12) 

where (-, ■) is the usual duality product for distributions on the sphere S. 

Since it is sometimes more convenient to work with mean zero functions (in order 
to use the main lemma of the previous subsection), we reformulate this problem in 
another framework. We set / = 1 + (7 so that / is a weak solution if and only 
if g e L\{0,T),H'+\S)) H\{0,T), H'-\E>)) with, for almost all t G [0,r], and 
for all h G H-'+^{S), 

{dtg, h) = -a{V^g, V^h) + (n - l)J[g] ■ J[h] + {g, J[g] ■ V^h). (13) 

That makes sense to look for a weak solution with prescribed initial condition 
in H^, since it always belongs to C([0,T], //'*(§)), as stated by the following propo- 
sition. 

Proposition 2.1. If g e L'^{{0,T), H'+\S)) n H\{0,T), H'-\S)) , then, up to 
redefining it on a set of measure zero, it belongs to C{[0,T], H^{§)) , and we have 

max \\u{t)\\l^ [ \\u\\l^^, + \\dtu\\l,_,, 
where the constant C depends only on T. 
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The proof in the case s = is the same as in j8], Thm 3, §5.9.2. To do the 
general case, we apply the result to {—A)'^g. 

Theorem 1. Given an initial probability measure /o in H'^(E>), there exists a unique 
weak solution f of such that /(O) = /q. This solution is global in time (the 
definition above is valid for any time T > 0). Moreover, f G C°°((0,+oo) x §), 
with f{t,uj) > for all positive t. 

We also have the following instantaneous regularity and uniform boundedness 
estimates (for m G N, the constant C depending only on a,m, s), for all t > 0: 

\\fmu^^c(^i + j^^ ii/oii^.. 

The proof consists in several steps, which we will treat as propositions. We first 
use a Galerkin method to prove existence on a small interval. We then show the 
continuity with respect to initial conditions on this interval (and so the uniqueness). 
Next, we prove the positivity of 1 + g for regular solutions. This gives us a better 
estimate of J[g]. Repeating the procedure on the following small interval, and so 
on, we can show that this extends to any t > 0. Regularizing the initial condition 
give then global existence in any case. 

We finally obtain the instantaneous regularity and boundary estimates by de- 
composing the solution between low and high modes. 

For the proof of all propositions, we will denote by Co,Ci,... some positive 
constants which depends only on s and a. We will also fix one parameter K > 
(which will be a bound on the norm of initial condition), and denote by Mq, Mi, . . . 
some positive constants which depends only on s and a, and K. 

Proposition 2.2. Existence: Galerkin method. 
We set 

where the constant Ci and C2 will be defined later. 

// ll^'ollifs ^ K, then we have existence of a weak solution on [0, T] satisfying (fT3l) . 
uniformly bounded m L^{{0,T), H'+\E,)) H H\{0,T), H'-\E>)) by a constant Mi. 

Proof. We denote by P/v the space spanned by the first (non-constant) eigen- 
vectors of the Laplace-Beltrami operator. This is a finite dimensional vector space, 
included in Hp{§>) for all p, and containing the functions of the form oj ^ V ■ u (see 
Appendix lA. 1 1 for more details). 

Let g'^ G C^{I,Pn) be the unique solution of the following Cauchy problem, 
defined on a maximal interval / C M_|_ ("non-linear" ODE on a finite dimensional 
space) : 

= U^iaA^g^ + (n - 1)(1 + g^) u ■ J[^/^] - J[^/^] ■ V.<?^), 
^7^(0) = n^((/o), 

where IlAr is the orthogonal projection on P/v- The first equation is equivalent to 
the fact that for any h G P/v, we have 

^ig'', h) = -a{V^g\ V.h) + {n - I) Jig''] ■ J[h] + {g'' , Jig''] ■ V^h). (15) 
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The goal is to prove that [0,T] G I and that there exists an extracted sequence Nk 
such that, as /c — >■ cx), 

• g^'' converges weakly in L^((0, T), if^"^^(S)) to a function g, 

• dtg^'^ converges weakly to dtg in L^((0,T), H^[E>)), 

• J[g^''] — 7- J[g] uniformly. 

We have that (— A)''^'^ G Pn, so we can take it for h, put it in (fTSl) and use the 
second part of Lemma [U to get: 

+ cr\\g''\\%,. ^ Co\J[gV\\ls + {n- l)V[^^]r (16) 
^C^yWl^il + C^yUs). (17) 

Indeed, any component of u belongs to any H^^, then J[g^] = {u,g^) is controlled 
by any if* norm of g'^ . 

Solving this inequality, we obtain for ^ t < Cf^ ln(l + (C2||n^(5fo)||ij3)~^)5 

.;v„ ^ lin^(^?o)b. 



\9 \\h^ 



e^^^^-C2\\U^{go)\\Hs{l-e~^^^)' 



Then we have ||g'^(t)||^3 ^ 2||(/o||^s for all t in [0,T]. There is no finite-time blow 
up in [0, T], then the ODE ffT^ has a solution on [0, T], for any N & N. 

Now we denote by AIq a bound for | on [0,T]. The inequality (fTUj) gives 



1 d_ 
2dt' 



/^lli + ^ll^?''llW^(l + ^o)C3|k||^.. 



Solving this inequality, we get for t G [0 , T] 

\\9 ll^s + o^ / \\9 \\h^+^ ^ WdoWfj,^ 
Jo 

We then use the ODE (ITSj) and this estimate to control the derivative of g. 

Taking h G L'^{{0, T), H~'^~^^{E>)), and integrating the equation in time, we get 

Jo 

Then we can take = K^e^^+^o^'^^^ max(a-\ C4 + Mq), and we get that g^ is 
bounded by Mi in L\{0, T), H'+\E)) f] H\{0, T), H'-\S)). 

Now, we just need estimates for ^J[g^]. We can take h = u-V in the ODE ( ITSl) 
and use the tools given in the beginning of this section. We get 

n — 1 



n 



(1 - an)J[g''] - / (Id - a; ® uj)J[g'']g'' du, 



^(C5 + MoC6)||^7o|lH^e^(i+*'»)^3T, 
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Indeed, again, since any component of Id — ui ^ u is in H we can control the 
term /g(Id — u ^ ^)g'^ da; by any H"^ norm of g'^ , uniformly in and in t G [0, T]. 

In summary if we suppose that go is in iJ*(S), for some s G M, we have that g^ 
is bounded in L\{0,T), H'+\E>))r) H\{0,T), H'~\S)), and that J[g^] and ^J[c/^] 
are uniformly bounded in and t G [0, T]. 

Then, using weak compactness and the Ascoli-Arzela theorem, we can find an 
increasing sequence A^^, a function g G L'^{{0,T), H'+^{S)) n H^{{0,T), H'-^{S)), 
and a continuous function J : [0,T] — > such that, as A; — oo, 

• J[g^''] converges uniformly to J on [0,T], 

• g^^ converges weakly to g in L'^{{0,T), H'+\E>)) and in H\{0,T), H'-^{§)). 

The limit g is also bounded by Mi in L^{{0,T), H'+\S)) D H\{0,T), H'-\S)). 
Then, since we have ip{t)uj{g^*' — (/) da; dt — )■ for any smooth function 93, 

we get Lp{t){J[g] — J) dt = and so J = J[g]. 

For a fixed h G Pm passing the weak limit in (fT^ (for A^^^ ^ M), we get for 
almost every t G [0,T] that 

V/i G Pm, (5^^?, /i) = -a{V^g, V^h) + (n - 1) J[(^] ■ J[h] + J[(/] ■ V^h). 

And this is valid for any M (except on a countable union of subsets of [0, T] of 
zero measure). By density (and using the first part of Lemma [I]), we have that g is 
a weak solution of our problem. 

Now for any h G H-'+\^), we have that {g^ {t) - YiN{.go):h) = J^{dtg^,h) 
is controlled by Mi-\/t||/i||^_s+i, uniformly in A^. So, passing the limit, we get 
that g{t) go in H-'+\E>) as t ^ 0. But since we know that g G C([0, T],H'{E>)), 
by uniqueness, we get g{0) = go- □ 

Proposition 2.3. Continuity with respect to the initial condition. 

Set T = ln(l + i_^2C2K ^' ^'^ (fT^ ■ Suppose we have two solutions g and g, 
with 11^(0)11^/^ K and ||/(0)||^. ^ K. 

Then there exists a constant M3 such that g — g is bounded in L^((0, T), i/*"'"^(§)) 
and m H\{Q,T),H^~\^)) by M3||^(0) - ^(0)||^.. 

This automatically gives uniqueness of a weak solution on (0, T) with initial 
condition go- 

Proof. Putting h = {—A)^g G in ffT^ . we do the same estimations as in the 

previous proposition. We have the same estimate as ( fT6|) -(fT71): 

IjtMls + a\\g\\l^,, ^ Co\J[g]\\\g\\l. + (n - 1)1%]^ (19) 

^CMls{l + C2\\gU^). (20) 

So if we set T = Cf^ln(l + (1 + 2C2K)~^), we can solve this inequality on [0,T], 
exactly as in fllSp . These solutions are then uniformly bounded in L^((0, T), if*+^(S)) 
and in H\{0,T), H'^\E>)) (by the constant Mi). 
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Taking u = g — g, and using ( fT3|) gives an equation for u: for almost all t G [0, T], 
for all h e H-'i^), 

{dtu, h) = -a{V^u, V^h) + {n- l)J[u] ■ J[h] + (n, J[g] ■ V^h) + {g, J[u] ■ V^h) . (21) 
Now we take h = {—AYu and use the first and second parts of Lemma [H to get 

l^^lMh + ^M\%s+i ^ (1 + Mi)Cs\\u\\l^ + Cr\\u\\j,4g\\Hs+A\i-^yu\\H~^ 

^M,{l + \\g\\^.+,)\\u\\l,. (22) 

Gronwall's lemma gives then the following estimate: 



\u 



h\\%s+i ^ \\uo\\%, exp (^M2 (1 + ||^||^,+i)^ 



Using dn}, we get that u is bounded in L'^{{0,T), H'+\S)) n H\{0,T), H'-\S)) 
by a constant M3 times ||m(0)||^s. □ 

Proposition 2.4. Positivity for regular solutions (maximum principle). Suppose 
that go is in i7'*(§), with s sufficiently large (according to the Sobolev emheddings, 
so s > is enough) so that the (unique) solution belongs to C^{[0,T],C^(S)). 
Here T is defined as in (HM . with K = WgoWn'^- 9'^ ^'^^/c to the original formula- 
tion f = 1 + g. Then f is a classical solution of ([3]). 

// /o is nonnegative, then f is positive for any positive time, and more precisely 
we have the following estimates, for all t G (0,T] and u ^ S (if /o is not equal to 
the constant function 1 ): 

e-("-i)/o inin/o < f{t,uj) < e^^-^^^o* l-^I/H max/o. (23) 

Proof. Since the solution is in C''([0, T], C^(S)), we can do the reverse integration by 
parts in the weak formulation (fT2|) . We get that, as an element of L^((0, T), H^~^{E>)), 
the function dtf is equal (almost everywhere) to aA^^f — ■ ((Id — w (g> oj)J[f]f), 
which is an element of C'^([0,T] x S). So up to redefining it on a set of measure 
zero, the function / belongs to ^^([O, T], C(§)) H C°([0, T],C^{S)), and satisfies the 
partial differential equation. 

Applying the chain rule and using the tools given in the beginning of this section, 
we get another formulation of the PDE ([S]): 

dj = aA^f - J[f] ■ V^/ +{n- l)J[f] ■ u f. (24) 

The next part of the proposition is just a classical strong maximum principle. 
We only prove here the left part of the inequality, the other part is very similar, once 
we have that / is positive. 

Suppose first that /o is positive. We denote by T > the first time such that 
the minimum on the unit sphere of / is zero (or T = T if / is always positive). 
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Then we have for t G [0, T], that dtf ^ aA^f - J[f] ■ V^f - {n - 1)| J[/]|/. If 
we write / = f e~^'^~^^ l-^I/H, we get 

dJ^aAj-J[f]-Vj. (25) 

Then the weak maximum principle (see [8], Thm 8, §7.1.4, which is also valid on the 
sphere) gives us that the minimum of / on [0,T] x S is reached on {0} x S. That 
means that we have a non-strict version of the left part of the inequality (l23|l : 

Vt G [0,f],Vcj G ^J{uj,t) ^ e^(""^)^ol^[/]linin/o. (26) 

Consequently, we have that min§ /(T) > and so T = T. If now /q is only 
nonnegative, take /q = and by continuity with respect to initial condition, 

inequality ( l26l) is still valid. That gives that / is nonnegative on [0,r], and conse- 
quently we have that inequality f l25l) is valid on [0,T]. 

Now we can use the strong maximum principle (see [S], Thm 11, §7.1.4), which 
gives that if the inequality ( l26l) is an equality for some t > and a; G S, then / is 
constant on [0,t] x S. So /q is the constant function 1. □ 

Proposition 2.5. Global existence, positivity. Suppose Jq is a probability measure 
belonging to H'^(S) (this is always the case for s < according to Sobolev embed- 

dings). Then there exists a global weak solution of which remains a probability 
measure for any time. 

We remark that the uniqueness of the solution on any time interval remains by 
Proposition 12.31 

Proof. We first prove this proposition in the case s > 

We define a solution by constructing it on a sequence of intervals. 

We set Ti = ln(l + i4,2C2||goH f|T^ . This gives existence to a solution g 

in C([0, Ti], i7''(S)). By induction we define Tk+i = Tk + (^ln(l + i+2C2||g(rfc)||jjJ ' 
which gives existence to a solution g G C{[Tk, T^+i], if*(S)). 

So we have a solution on [0, T], provided that T ^ Tk for some integer k. 

Now by the previous proposition, this solution / = 1 + g is nonnegative. We 
obviously have \ J[g]\ = \J[f]\ ^ Is\^\f — 1- Then we can do better estimates, 
starting from (fT^ : 

lft^\9\\ls + cr\\g\\l.,, ^ Co\J[9]\\\g\\l. + (^- ^)Vl9]\' 

^CsMl^. (27) 

Then, Gronwall's lemma gives us that ||(7(Tfc)||^3 ^ ll5'o||i/se'"*'^'=. Suppose now that 
the sequence (T^) is bounded, then ||(7(Tfc)||^3 is also bounded. By the definition 
of Tfc_|_i, the difference T^+i — T^ does not tend to zero, which implies that the 
increasing sequence (T^) is unbounded, and this is a contradiction. So we have 

that Tk oo, and the solution is global in time. 
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Now we do the general case for any s. Take a sequence of elements of H^~^^ 
converging to qq in H'^, and such that /q = 1 + are positive functions. Let be 
the solutions associated to these initial conditions. 

Then we have the same estimates as before, since we still have \J[g]\ ^ 1, solv- 
ing (177|) gives 



Jo 



So we can now study the difference u = g^ — g\ as in (l^T]) . fl2^ . which satisfies, for 
any h e i^-"(S), 

{dtu,h) = -a{V^u,V^h) + {n-l)J[u]-J[h] + {u,J[g'']-V^h) + {g^,J[u]-V^h). (28) 
We take h = (— A)*ii and use the first and second part of Lemma [T] to get 

1 d 
2dt' 

^Cw{1 + \\9\\hs+i)\\u\\1,. (29) 
If we fix T > 0, Gronwall's lemma gives then the following estimate: 



+ f^lhllis+i ^ C'glliillJ^, + C7||w||^.||5i''||^.+i||(-A)'*m||^_. 



Since ||g'o||ii-s is bounded (because g^ converges in H^), together with f l2T]) . we finally 
get that u is bounded in L2((0, T),H'+\E>)) nH\{0,T), H'-\S)) by a constant Ct 
times ||m(0)||j^3. This gives that g^ is a Cauchy sequence in that space, and then 
it converges to a function g, which is a weak solution of our problem (by Proposi- 
tion [2]T1 we have that g{0) = go). This is valid for any T > 0, so this solution is 
global. 

If we take ip in C°°(S), since f^{t) = 1 + g^{t) is a positive function with mean 1, 
we have that 

-IIV'lloo = {f{t),-M^) ^ {f\t),^) ^ (fit), llv^lloo) = Moo. 

Then passing the limit gives \{g{t),if)\ ^ ||v^||oo- Furthermore we have {f^{t)^ 1) = 1 
so (/(t),l) = 1, and if is a nonnegative function, then {f''{t),if) ^ and we 
get {f{t),if) ^ 0. This gives that f(t) is a positive radon measure with mass 1, 
which is a probability measure. □ 

Proposition 2.6. Instantaneous regularity and boundedness estimates. If /o is a 
probability measure, then the solution f belongs to C°°((0,+oo) x §), is positive for 
any time t > 0, and we have the following estimates, for a// s G M and m ^ 0: 



t 

where the constant C depends only on a, s, and m. 

In particular we have that for to > 0, f is uniformly bounded on [to,+oo) in 
any norm. 
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Proof. Suppose /o G H'{S), and fix t > 0. The solution / is in C([0, +00), 
and in L^{lo,t), H'+\S)). Then there exists s < t such that f{s) G H'+^{E>). So 
we can construct a solution belonging to C([s, +00), But this solution 

is also a weak solution in L^{{s,T), H'+\S)) n H~\{s,T), H'-\S)), for all T > s 
so by uniqueness it is equal to /. Then / belongs to C([t, +00), if*"'"^(S)). Since 
this is true for all t > 0, then / belongs to C{{0, +00), H^~^^ {§>)). We can repeat 
this argument and have that / belongs to C{{0, +00), H^{S)) for any p, and is a 
positive classical solution, by Proposition 12.41 Using the equation, differentiating in 
time gives that it is also in C'^((0, +00), if^(S)) for any p and any k, so, by Sobolev 
embeddings, it is a C°° function of (0, +00) x §. 

Since we have positivity, we can have estimates for any of the modes of f = 1 + g. 
Let us denote the orthogonal projection of / on the first eigenspaces of the 
Laplacian, and = f — the projection on the other ones (high modes). 

We have a Poincare inequality on this space: llfi'^ll^^ ^ (Af+i)(jv+n-i) ll^^^ll jj-s+i 
(we recall that the eigenvalues of —A are given by £(£ + n — 2) for £ G N). We use 
the estimate fll9p : 

\^Mh^ + ^ll^lli-^ ^ Co\J[g]\\\g\\\. + {n- 



^ xwrmr^Mls^^ + - + C4f - (so) 

Now we have, since / is a probability measure, that 

lir - = / (-A)V^/du; ^ ||(-A)7^1|loo ^ K^lir - 

the last inequality being the equivalence between norms in finite dimension. Dividing 
by this last norm, this gives that the low modes of / are uniformly bounded in time 
by a constant K^. Then we have, taking sufficiently large, 

1 d 2 II ||2 ^ n 

Now multiplying by t this formula at order s + 1, we get 



and finally 

Together with Poincare inequality, solving this inequality gives us 

Mis + ^ ll^?o|li.e-("-^)t* + C^s{l + t). 

So we have the result for = 1 + llfi^ll^^, and m = 1: 

\\m\\w^c(^i + ^^ \\M\%. 

Then we apply this inequality between and |, and the inequality at order m 
between | and t to get the result at order m+1. The case where m is any nonnegative 
real also works, by interpolation. □ 
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This last proposition ends the proof of Theorem [H Let us do here two small 
comments concerning the analyticity of the solution and the limit case with no 
noise: cr = 0. 

Remark 2.1. Analyticity of the solution. We can show, as claimed in f^, ^ that 
at any time t > the solution is analytic in the space variable. The idea is to 
show, following (based on J^, fW^). that the solution is in some Gevrey class of 
functions, defined by a parameter depending on time. This class is a subset of the set 
of real analytic functions on the sphere. More details and a complete proof are given 
in Appendix ] A. ^ We could have directly dealt with this classes of functions instead 
of working in the Sobolev spaces, but we will not need these properties of analyticity 
in the following. In any case, to prove analyticity we need the initial condition to be 
in H~^~{E>), so this study of instantaneous regularization was necessary. 

Remark 2.2. Case where a = 0: no noise. The proof is also valid, except that 
the solution belongs to L°°((0, T), n H^{{0,T), H'-\E>)) if the initial condi- 

tion is in H'^{E>). By an optimal regularity argument, we can get that a solution is 
in fact in C([0, T], if^(S)). The nonnegativity argument is then also valid, and so 
the solution is global. Obviously, we do not have the instantaneous regularity and 
boundedness estimates. 

3 Using the free energy 

In this section, we derive the Onsager free energy ([6]) for Doi equation ([3]), and use 
it to get general results on the steady states. 

3.1 Free energy and steady states 

We rewrite the equation ([3]): 

dtf = Qif) = V. ■ (aV./ - V^co ■ J[f])f) = V. • (/V.(aln/ - u ■ J[f])). 

Since any solution is in C°°((0,+oo) x S), and positive for any t > 0, there is 
no problem with using In /, and doing any integration by parts. We multiply the 
equation by crln / — a; ■ J[f] and integrate by parts, we get 

/ dJ{a\nf-u-J[f])du = - [ f\VUcr\nf-u-J[f])\'du. 

Since the left part can be recast as a time derivative, this is a conservation relation. 
We define the free energy J^{f) and the dissipation term by 

:F{f) = a [ /ln/-i|J[/]p, (31) 
V{f)= [ /|V.(aln/-a;- J[/])|^ (32) 
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and we have the following conservation relation: 




(33) 



We define a steady state as a (weak) solution which does not depend on time. 
Here are some characterizations of the steady states. 

Proposition 3.1. Steady states. The steady states of Doi equation ([3]) are the 

probability measures / on S which satisfy one of the following equivalent conditions. 

1. Equilibrium: f G ^^(S) and Q{f) = 

2. No dissipation: f e C\S) and V{f) = 

3. The probability density f G C'^(S) is positive and a critical point of T (under 
the constraint of mean 1 ). 

4. There exists C G M such that a In / — J[f] ■ w = C. 

Proof. By definition, a steady state / is a solution independent of t. Since it is a 
solution, it is positive and C°°, and we get that Q{f) = 0. By the conservation 
relation f l33l) . we get that ^J-" = 0, so P(/) = 0. Since it is positive, we get 
that Inf — u- J[f]) = 0, so there exists C G M such that cr In / — J[f] ■ uj = C. 

Now we do a variational study of J-" around /. We take a small perturbation f + h 
of / which remains a probability density function (which means that J^h = 0). 

We can expand the function x ^ xlnx around /, since f ^ e > 0, and we have 



which means that / is a critical point of J-". So / satisfies the four conditions. 

Conversely if / G C^(S) and Q{f) = 0, then / is obviously a steady-state. 

If (Tin / - J[f] ■u = C, then / G C\E>) and Q{f) = 0. We will show that the 
second and third conditions reduce to this fourth condition. 

Doing the above computation around a positive / G C^{S) gives that if / is a 
critical point for the free energy, then j^h{alnf — J[f] ■ u) is zero for any h with 
mean zero. This is exactly saying that cr In / — J[f] ■ u is constant. 

Finally if we suppose / G C^(S) and 'P(/) = 0, at any point G § such 
that /(wq) > we have that V(o"ln/ — J[f] ■ cj) = on a neighborhood of Uq. 

The function ip defined by (p{uj) = cr In / — J[f] ■ u is then locally constant at any 
point where it is finite, so Lp~^{{C}) is open in § for any C G M. 

Now if (p{cdk) = C*, with Uk converging to u^o, then f{ujk) = exp(^^i^^^^l^). 
Passing to the limit, we get that f{u^) = exp( '""^'^[^^''^°° ), which gives (p{uj^) = C. 
So ip~^{{C}) is closed. 

Since / is not identically zero, there exists C G M such that Lp~^{{C}) 7^ 0, and 
by connectedness of the sphere, we get if^^{{C}) = §, so aln/ — J[f] ■ oj = C. □ 



Hf + h) = a f{f\nf + h\nf + h)- \\J\ff - J[f] ■ [uh + Oml) 
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3.2 LaSalle principle 

We give here an adaptation of LaSalle's invariance principle to our PDE framework. 

Proposition 3.2. LaSalle's invariance principle. Let fo be a probability measure 
on the sphere §. We denote by J^oo the limit of J^{f{t)) as t ^ oo, where f is the 
solution to Doi equation ([3]) with initial condition /q. 

Then the set S^o = {/ ^ C°°(§) s.t. V^f) = and J^{f) = J^oo} is not empty. 

Furthermore f{t) converges in any norm to this set of equilibria (in the fol- 
lowing sense): 

lim inf = 0. 

Proof. First of all J^{f{t)) is decreasing in time, and bounded below by — |, so J^^o 
is well defined. 

Let {tn) be an unbounded increasing sequence, and suppose that f{tn) converges 
in H^{Ei) to /oo for some s G M. We first remark that /(t„,) is uniformly bounded 
in H^~^^^{E>) (using Theorem [1]), and then by a simple interpolation estimate we get 
that ||/(t„) - /(t™)||^.^, ^ Wfitn) - f{tm.)\\Hs\\f{tn) - /(^„^) || and /(t„) also 
converges in H^~^^{E>). So /oo is in any H^{E>). 

We want to prove that P(/oo) = 0. Supposing this is not the case, we write 

V{f) = a' [ + J[f] ■ /(Id - a; ® uj)f J[f] - 2aJ[f] ■ [ V./ 

Js J Js Js 



/ 



a' I + (1 - 2(n - l)a)\J[f]\' - [ {uj ■ J[f])'f. (34) 

Js 



Now we take s sufficiently large such that //*(§) C Loo(S>) r\H^{E>). If /oo is positive, 
then V, as a function from the nonnegative elements of if*(S) to [0, +oo], is contin- 
uous at the point /oo. In particular since D(/oo) > 0, there exist 6 > and M > 
such that if ||/ — fooWn^ ^ ^, then we have X'(/) ^ M. We want to show the same 
result in the case where /oo is only nonnegative. We define 

Vsif) =a' I + (1 - 2{n - l)a)\J[f]\' - f {co ■ J[/])V. 

Js J ~r ^ Js 

We have that by monotone convergence that ^^^(/oo) converges to D(/oo) as e — )■ 0. 
So there exists e > such that T>^{foo) > 0. Now by continuity of at the 
point /oo, we get that there exists 6 > and M > such that if ||/ — /oo||h= ^ 
then T>^{f) ^ M. And the fact that !'(/) ^ T^eif) gives the same result as before. 

Now since dtf is uniformly bounded in H'^ (for t ^ ti > 0), there exists r > 
such that if \t-t'\ ^ r, then \\f{t) - f{t')\\Hs ^ f- We take then N sufficiently large 
such that II /(t„) - /oolk- ^ I for all n ^ iV. 

Then we have that for n ^ N, 'D{f) ^ M on [tni^n + t]- Up to extracting, we 
can assume that t„_|_i ^ t„ + r, so we have 

nfM) - Hf{tN+p)) = / V{f) > prM. 

JtN 
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Since the left term is bounded by J-'(/(tAr)) — J^oo, taking p sufficiently large gives 
the contradiction. 

Now if we suppose that for a given s the distance (in H'^ norm) between f{t) 
and £00 does not tend to 0, we get e > and a sequence t„ such that for all g G £^00; 
we have ||/(tn) ~ fl'||//= ^ ^- Since is bounded in H'^^^iS), by a compact 

Sobolev embedding, up to extracting we can assume that is converging in H'^(S) 
to /oo- By the previous argument / G C°°{S) and we have P(/oo) = 0. Obviously 
since is decreasing in time we have that J^(/oo) = J^oo- So /oo belongs to £"00, 

and then — /ooWh" ^ £ foi' all n. This is a contradiction. 

Since the distance between f{t) and £"00 tends to 0, obviously this set is not 
empty. □ 



3.3 Computation of equilibria 

Define, for a unit vector G S, and n ^ the Fisher- Von Mises distribution with 
concentration parameter k, and orientation Q by 

M^f, a; = -f— . (35) 

Note that the denominator depends only on k. We have that the density of M^fj 
is 1, and the flux is 

J[M^n] = Y = c{K)n, (36) 

Jg exp(«:a; ■ izjdo; 

where 

£ cos ^ e-^ ^"'^ ^ sin""2 ^ 
'^^''^ " /;e--^sin"-2^d^ • ^^^^ 

If / is an equilibrium, cr In/ — J[f] -u is constant, and then / = C exp(cr"^ J[/] -u). 
Since / is a probability density function, we get / = M^q with = cr~^ J[/] (in the 
case where \J[f]\ = 0, then k = and we can take any Q, this is just the uniform 
distribution). Finally with (l36l) we get J[f] = c{k)Q, which gives the following 
compatibility condition 

c{k) = an. (38) 
We give the solutions of this equation in a proposition. 

Proposition 3.3. Compatibility condition 

• If a ^ ^, there is only one solution to the compatibility condition: k = 0. The 
only equilibrium is the constant function / = 1. 

• If a < ^, the compatibility condition has exactly two solutions: k = and one 
unique positive solution, that we will denote K{a). Apart from the constant 
function / = 1 (the case k = 0), the equilibria form a manifold of dimen- 
sion n — 1: the functions of the form f = Mi^(^„-jq, where Q ^ Ei is an arbitrary 
unit vector. 
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Proof. Let us denote cr(/«) = A simple Taylor expansion gives (t(k) — S-^ K Since 

the function a tends to as k — )■ +00 (because c{k) ^ 1), it is sufficient to prove that 
it is decreasing. Indeed the function is then a one-to-one correspondence from 
to (0, and the compatibility condition for A€ > is exactly solving a = a^n). 

But we have (after one integration by parts) that ?'(«:) = -^(1 — na{K) — c(k)^), 
which, by the following lemma is negative for k > 0. □ 

Lemma 2. Define /3 = c(/t)^ + na{K) — 1. Then for any k, > 0, we have /3 > 0. 

Proof Define [7(cos^)]^ = /^^ 7(cos^) e''"°'^^ sin""^ ^d^. 

rp, , u A R K[cose]l + n[cos9l[ll-K[l]l 

ihen we have by definition p = -— . bo we only 

have to show that the numerator is positive. We will prove in fact that the Taylor 
expansion of this term in k has only positive terms. 

We have, if we denote = -^2^ cos^^ 6 sin"~^ 9d9 ^ 0, 

00 00 
= XI ^p'^''' and [cos 6*], = J](2p + 2) ap+iK^P+^ 

p=0 p=0 

Now doing an integration by part in the definition of a^+i, we get 

_ 2p + 1 / ap 
" V(2p+l)(2p + 2) " ' ' 

which gives 

We have, for k > 0, 



= + 2) ap+,{2q + 2) a,+i + J] n(2p + 2) a^+i 

k=0 \p-\-q=k—l p-\-q=k 



Qjq QipQjq | Ki 



fc=0 \p+<j'=fc,p^l p+q=k 



2k+l 



2p+n I "'P"-1 



2(j+n 2p4 

fc=0 \p+g=fc 
00 / 



2fc+l 



fc=0 \p+q=k 

00 



= E E HSKifcj''.s h"" 

fc=0 \p+q=k / 

So we finally get 

= ( X] '^P*^^ ] X] ( (2p+n)(2g+n) ^P*^? ) 

\p=0 / fc=0 \p+ij=fc / 

which gives that /3 > when k > 0. □ 
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Remark 3.1. We can do another proof, following an argument of [2^ . which does 
not need to compute explicitly p. 

The idea is that we compute a" = (n — 1)^ — 2a{a — /3), so we see (except in 
the case n = 0) that if a' = — ^ = 0, then a" < (indeed, we will easily see in 
that a — f3 is positive). For the case k = 0, we can compute the Taylor expansion 
of a up to order 2; = ^ — n^(n+2) ^^ 0{k^). So we have that any critical point 
of a is a maximum. Since there is a local maximum at n = then the function is 
decreasing. 

We can have an asymptotic expansion of the order parameter c{k{cf)) as a reaches 
the critical value ^. Indeed we have that c" — ;^ ~ ~ ra2(n+2) ^('^)^ expansion 
of a in the previous remark. So 

c(K(a)) ~ ^K{a) ~ y/(n + 2)(i-a) as a ^ 1. (40) 

Proposition 3.4. Minimum of the free energy 

• If a ^ ^, the minimum of the free energy is 0, only reached by the uniform 
distribution. Any solution converges to the uniform distribution in any 
norm. 

• If a < j^, the minimum of the free energy is negative, only reached by any 
non-isotropic equilibrium M^i^^-^^i. 

Proof. By LaSalle principle (Proposition 13. 2p . we have that 

min J^if) = niin J^if)- 

/ec°°{§),/>0 /ec=°(§),/>o,©(/)=o 

Indeed for any positive initial condition / in C°°{S), there exists an equilibrium f^o 
such that F{foo) = ^oo ^ J^if)- This gives 

inf -F(/)= inf J^{f). 

/ec°°(§),/>o /6C°°(s),/>o,©(/)=o 

Since the set of equilibria is compact (either a single point or one point and a 
manifold homeomorphic to S), this infimum is a minimum. 

Furthermore, if /o is not an equilibrium, then ^^(/o) > 0, and then J^{f{t)) is 
decreasing in the neighborhood of t = 0. So the minimum of J-" cannot be reached 
for /o. 

In the case a ^ this gives the result since the only equilibrium is the constant 
function 1. By LaSalle principle, we also get that the solution is converging to in 
any H'^ norm. 

In the case a < we have that ^-"(1 + eoj ■ Q) ~ ^{a — ^)e'^ for a fixed unit 
vector f2 G S, so there exists /o such that J^{fo) < 0. Then the uniform distribution 
cannot be a global minimizer. Since J-'(MK(o-)n) is independent of Q, we get that this 
value is the minimum. □ 
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4 Convergence to equilibrium 



In this section, we establish and study the convergence of the solution to an equi- 
librium for any initial condition, in the three different regimes, depending wether a 
is greater, less, or equal to 

4.1 A new entropy, application to the subcritical case a > K 

In this section we derive a convex entropy, which shows global decay to the uniform 

distribution in the case a > -. 

1 ~ 
We define on H~^~{E>) the norm || -H^.n^ by = J^d^n-id^ where the 

conformal Laplacian A„_i is defined by ([7]). This norm is equivalent to || • ||^_n^. 

We also define || ■ ||^_n^ by ||(7||l_„-3 = J^AgA~^-^^g, and this norm is equivalent 

to the II ■ il . n-3 norm. 

Taking h = A~\g in the weak formulation ( fT3|) . and using the last part of 
Lemma [H we obtain a conservation relation: 

~\\g\\l = -a\\g\\l „-3 + , ^ ,J J[9]\'. (41) 
2dt" "h-— ii^iI//-^ (n-2)!' ^ ' 

We remark that this is a conservation law between quadratic quantities, as it would 
be the case for a linear equation. 

Since the component of g on the space of spherical harmonics of degree 1 is 
given by nu ■ J[g], a simple computation shows that the contribution to Hf^Hl „-i 

H 2 

of this component is equal to ^^^-ly. I"^!^]!^- Then the last term of the conservation 
relation f l4ip is bounded by -^^H^fUl ,,.! . Together with the Poincare inequal- 
ity ll5'll~_i-3 ^ ~ 1)||5'|| ~_n-3 ) we get the following estimate: 

This gives in the case cr > ^ an exponential decay of rate {n — l)(cr — ^) for the 



gWfr-^-^ ^ IkoL-ii^ exp(-(n - l)(a- ^)t). 



In the general case, if fo e H'{^) with s > we use the estimate ([30]): 

Now we have, since / is a probability measure, 

the first inequality being the equivalence between norms in finite dimension. For 
any £ < ^, taking sufficiently large, together with Poincare inequality we get 

liM% + - - ^)M% < C||^o||^.:^e-("-)(--^^ 
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where the constant C depends only on s. 
Solving this equation, we get 

< llnJI^ -2(n-l){a-e)t . C 11 ||2 2(n-l)(<7- i)t 

Taking for example £ = since s > — we get 
In summary, we have the following theorem: 

Theorem 2. New entropy. For a given probability density function f , we define the 
quantities nif) = \\f-l\\l_^^, and Vif) = 2a\\f - -^\J[f]\^. 

We have a conservation relation, for any solution f of Doi equation 

^n{f)+V{f) = 0. (42) 

When a ^ ^, the term T>[f) is nonnegative, so the new entropy H{f) is decreasing 
in time. 

Furthermore, if a > ^, then in any Sobolev space H^(E>) with s ^ ^'^^j '"^c have 
global exponential decay of the solution to the uniform distribution, with rate given 
by (n-l)(a-i). 

More precisely there is a constant C depending only on s such that for all initial 
condition /o G iif''(S), we have 

Let us do a small remark here. Actually this conservation relation is true for 
any solution, without any positivity condition. We only need the mean of / to be 1. 
And since we have existence and uniqueness in small time for any initial condition, 
with the same instantaneous regularity results (only valid for a short time existence), 
we get that the solution belongs to i/~^(S) at some time. But the conservation 
relation gives then that we have a global solution. So we can state a stronger theorem 
of existence and uniqueness: 

Theorem 3. Given an initial condition /q in H^(E>) (not necessarily nonnegative), 
there exists a unique weak solution f of such that /(O) = /q. This solution is 
global in time (the definition \2.1\ is valid for any time T > 0). Moreover, f is a 
classical solution, belonging to C°^((0,+oo) x §) (and even analytic in space, see 
Appendix ] A. 

Remark 4.1. In this case, we do not have any uniform bound on H''^{E>), and we 
can derive the same existence theorem for the case cr = (see Remark \2.^) . but only 
for the case s ^ (which does not include all radon signed-measures) . 

Another remark is that if we change the sign in front of the alignment term in Doi 
equation ([3]) (taking K{u,u}) = u-ui, every particle tends to go away from the mean 
direction), then we can derive a conservation relation in the same way. But here 
the "dissipation term" is V{f) = 2a\\f - + ^^\J[f]\^ > 2a(n - l)n{f), 

without any condition on a > 0. So in any Sobolev space H^{S), with s > — 
we have global exponential decay of the solution to the uniform distribution, with 
rate in — l)cr. 
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4.2 Study of the supercritical case < ^ 

In this section, we fix a < ^ and we study the behavior of a solution as t — )■ +00. 
We will write k for K{a) and c for c(k(o")). We first establish that the limit set of 
equilibria £^00 given by LaSalle principle (Proposition 13. 2p depends only on the fact 
that J[/o] is zero or not. 

Proposition 4.1. // J[/o] = then £^ is reduced to the uniform distribution. 
Equation ([3]) becomes the heat equation. We have exponential decay to the uniform 
distribution with rate 2na in any H^{E>). 

If J[/o] ^ then J[f{t)] ^ for all t > 0. The limit set S^o = {M^n,n E §} 
consists in all the non-isotropic equilibria. Furthermore, we have for any s G M, 

lim\\f{t)-M,nit)\\Hs=0, (43) 
where f2(t) = is the mean direction of f{t). 

Proof. First of all, we write the equation for multiplying equation ([3]) and 

integrating on the sphere. We get 

^ J[/] = -a{n - l)J\f] + Q(Id -u®uj)fdu}j J\f] 

(l_(„_l)^)Id- [co®uf]j[f], (44) 



which can be viewed as a first order linear ODE of the form ^ J[/] = M{t)J[f]. 
The matrix M is a smooth function of time, so we have a global unique solution. 
Consequently, if J[/(to)] = for to ^ 0, then we have J[f{t)] = 0, for all t ^ 0, and 
equation ([3]) reduces to the heat equation. The distribution / has no component 
on the first eigenspace of the Laplace-Beltrami operator, and the second eigenvalue 
is 2n, so we have exponential decay with rate 2na in any H'^ norm. 

Now we suppose that J[fo] 7^ 0, so by the previous argument we have J[f{t)] 7^ 
for all t ^ 0. There are two possibilities for the limiting set, either the uniform 
distribution, or the set {Mf^Q,Q G S} (by Proposition 13. 4[ they do not have the 
same level of free energy). 

In the first case, by LaSalle principle, f(t) converges to the uniform distribution. 
Then the matrix M{t) = (1 — (n — l)cr)Id — J^u f converges to (n — 1)(^ — cr)Id. 
Using the ODE for J[f], we get 

l^tWW = Af] ■ Mmf] ^ {{n - l)(i -a)- e)\J[ff, 

for t sufficiently large. Taking e sufficiently small, we get that \J[f] \ tends to infinity, 
which is a contradiction. 

So we have that Soo = {Mi^q, Q G S}. Now suppose that \\f{t) — M^fK^t^Wn" does 
not tend to 0. We take t^ tending to infinity such that \\f{tn) — Mf^Q(^tn) Wh" ^ s > 0. 
By our LaSalle principle, there exists f2„ G § such that ||/(tn) ~ M^^q^Wh^ — )■ 0. Up 
to extracting, we can suppose that Qn ^00 G S, so f{tn) — M^nac in H'^{S>). In 
particular we have that J[f{tn)] — J- c{k)Qoo, and then Q{tn) ^oo- Then M^nitn) 
converges to M^n^, giving that ||/(tn) — ^Kn(t„)||iT'' 0, which is a contradiction. 

□ 



23 



Now we focus on the case J[fo] 7^ 0. We define Q{t) as in the previous proposition, 
and we will expand the solution around Mi^^{ty We first show the convergence 
in to a given equilibrium, with exponential rate, assuming conditions on the 

initial data. 

Proposition 4.2. There exists an "asymptotic rate" roofer) > satisfying the fol- 
lowing property. 

Suppose that MKQ(t) is uniformly bounded on [to, +00) by a constant K , 

with s > ^^"-'^^ , Then for all r < r^oio'), there exists il^o ^ S and 5,C > 0, such 
that if II /(to) - M^Q^to)\\L^ ^ we have 

\\f{t) - M^nJlL^ < C\\f{to) - M^mtjL^e-^^'-''^. 

The constants 6 and C depend only on a, s, K, and r. Moreover, as a ^ ^, we 
have that r^{a) ^ 2(n - \){\ - a) + 0((i - a)t). 

Proof. We first introduce some notations. When there is no confusion, we just 
write Q for Q{t), and we will always assume t ^ to- We write cos^ = 00 ■ Q. We 
denote by {■)M^i^ the mean of a function against the probability measure M^q. 

We have the following identities (we recall, by Lemma [21 that /3 = + ncr — 1 is 
positive) : 

(cos^^)M,n = 1 - (n - 1)CT, 
((cos^ - c)^)M,n = I - {n - l)a - = a - P > 0. (45) 

We can write / = (1 + h)Mi^Q, then we have {h)M^Q = 0. Since Q is the direction 

of Af] = ((1 + h)uj)M^n^ we get that {hu)M^a = {hco^9)M^n^- 

So we can do an expansion of the free energy and its dissipation in terms of h. 
Since we know that M^Q(^t) is a critical point of J-", we already know that the expan- 
sion of J-'((l + h)Mi^Q) — ^'{Mf.n) will contain no term of order and 1 in h. We 
get, using (EI]), 

^((1 + h)M^n) - HM.n) = a|(/i^)M.. - ||(Ma/.oI' + Oml)- 

Using Sobolev embedding and interpolation, we have (writing C for a generic con- 
stant, depending only on a, s, and K) 

||/-M«f,|U ^C||/-M,n||^-i ^C||/-M,nL2 K—. 

So since 1 — > | and / — M^f^ = HM^q, with M^q uniformly bounded below 
and above, we get that ||/i||^ = o{{h?) j^i^^) (and more precisely, for any e > there 
exists ?7 > depending only on e, cr, s, and K such that ||/i||i^ ^ ^(^^)A/„n as soon 
as (/i^)M«n ^ r/). We get 

Hf) - HM^n) = lW{h')M.,, - {hcose)lJ + o((/i2)m.J. (46) 
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We use the definition §^ oiV{f): 

D{f) = ((1 + h)\V{a\n{M,n{l + h)) - ((1 + h)u)M^n " ^)\'')m^ 



n 

((1 + /i)|V(aln(l + h)- {hcose)M^,, cos^)Hm.o 



^ (1 - ||/i||oo)(|V(aln(l + /i) - {hcose)M^a cos^)nM, 
Now we can derive a Poincare inequality of tfie form 

{\^9?)M.n^^n{{g-{g)M^a?)M^a- 

Indeed, we use tlie fact tliat M^q is positive and bounded: 

> min M^Q{n - 1) / - /§ g) 



(47) 



(45 



^(n-l)e-^'^((^7-((7)M.J')A./.o- (49) 

Actually this is a rough estimate, we have here ^ (n — l)e~^'', a more precise 
study of Ak could be done using separation of variable. The problem then reduces to 
finding the smallest eigenvalue of a one-dimensional Sturm-Liouville problem, but 
even in that case, we did not manage to find a better estimate for now. 
So we finally get 

V{f) ^ (1 - ||/i|U)A«([aln(l + h)- a(ln(l + /i))a/.,- {h cos 9) {cos 6 - c)f)M^n 
^ (1 - \\h\\^)A^{[ah - {hcose)M^,,icose- c) + Oi\\hfJ]^)M^^ 
^ (1 - \\h\\^)A^ia'{h')M^, - (/3 + a){hcose)l,J + OmD- 

With the same argument as before, we get that 

V{f) ^ A^{a'{h')M^,, - (/3 + a){hcose)l,J + o((/.2)m.J. (50) 

The goal is now to express the bounds in ( l50l) and ( H6l) as the sum of positive 
terms. Indeed, we expect to have a Gronwall's inequality which will give a rate of 
convergence. 



We set a = ■^^{hcos6)M^Q, and we write h = a{cos9 — c) + g. Using (j45 



we 



have that a is well defined since a — (3 > and we get {g)Ai^<^ = {gu)M^n ~ 0- 
Plugging {h^)M^n = (a - /3)a2 + {g^)M^n into (|46D and (EU]) gives 

^(/) - HM.^) = - fi)o? + a(/)M„J + o((/i2)m.J, (51) 

^ A,/3(/3(a - /3)a' + a(/)A/,J + o((/i2),,^J. 
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So for all r < A^/?, if {h?) m^q is sufficiently small, we have ^ r(J-'(/)— J-'(Mkq)). 

Using the conservation relation ( l42l) . there exists 5o > (depending only on a, s, K 
and r) such that if \\f{t) — Mi^Q(^t)\\L2 ^ 60, we have 

^[^(/) - T{M,n)] = -V{f) ^ -2r[^(/) - ^(M.^)]. 
Then we obtain, for all T, such that ||/ — M^^qW^-i ^ 5o on [to, T], 

and then, using the estimate ( 15T]) . we get that for t G [to,^], 

11/ - M^nh^ ^ Coif (to) - M«f,(i„)||L2e-''(*-*°). (52) 

So if we take 6 < ^ 60, and we start with ||/(to) — ^Kn(to)||L2 ^ 5, we get 
that 11/ — Mkq||l2 ^ ^0 on [^o,^] for all T ^ to- Otherwise, the largest of such a T 
would satisfy 60 = ||/(T) - M^f^(r)i|L2 ^ C5e-^(^-*«) < ^q. So the inequality ([52]) 
holds for all t G [to, +00). 

It remains to prove that Q{t) converges to some Q^o, if we want to have strong 
convergence to a given steady state. This is possible using the ODE satisfied by Q. 

Indeed, we have J[f] = cfl + {huj)M^n = (c+ a(o^ — Z^))^, and then 

^ J[/] = (c + a{a - /?)) + (a - Aa. 

So applying Id — (8> to the ODE (jH]) gives an ODE for f2, in terms of a and g. 
We get 



(Id - ^] ® f^)^ J[/] = -(Id - ® 1]) ( /" 
dt Vis 



w (g) a; / dw J[/] 



= -{c + a{a - (5)){\d-VL®VL)[{hco^eu)M^^ + (cos^a;)M,n]- 
Since ((cos 6' — c) cos6'a;)M«n {cos 9 u)m^q are parallel to f2, we get that 

(c + a(a-/?))^ = -{c + a{a- /3)) (Id- n(S)n){g cose u)M,n- 

Since (c + a(cr — is the norm of J[f], it is never zero, and we get (the notation C 
standing for a generic constant depending only on r, s, a and K) 



dn 



dt 



So we have exponential decay of ^ with rate r, in particular ^2 is converging to 
some floo £ S. More precisely. 



/oo 
\f\dt^C\\fito)-ALnito)\\L'^e- 



-rit-to) 



26 



Now we have that ||M«Q(t) - M«n^||i2 ^ C\Q{t) - Q^d (the function Q e'"^'^ 
from S to M is globally Lipschitz with a constant independent of a; G S). So we get 
the final estimation: 

\\f-M^nJ\L'^ ^ \\f-M,n\\L'^ + \\M,nit)-M^nJ\L^ ^ C\\fito) - M^nito)\\L^e-''^''''^ . 

So the proposition is true with roo(cr) = A^/? > 0. By the estimate ( 149|) . we know 
that Aft; ^ (n — 1)6^^". And by the expansions of c and k as cr — ^ given in (140|) . 
we get that roo(a) ^ 2{n - - cr) + 0((i - a)l). □ 



By Proposition 14.11 we have that /(t) — M^^Q(^fj tends to zero in any if*(S). So 
the hypotheses of Proposition 14. 2^ for any r < roo(o"), are satisfied for some to > 0. 

Once more, by interpolation and uniform boundedness on [to,+oo) of the 
norm, we have 

11/ - M^nJlHs ^ C\\f - M^nJlPWf - M^nJlL 

so taking p sufficiently large, we also get exponential convergence for the if^ norm, 
with rate r(l — 6) for any 6 > 0. 

Finally we have that for all r < roo{a) and s, there exists some time to and C > 
such that 11/ — M^q^Whs ^ Ce"''* for t ^ Iq. We can even get rid of the constant C 
since for any r < r and t sufficiently large Ce"''* ^ e"''*. 

4.3 Study of the critical case a = - 

For any a G (0, +oo) \ {-}, we have exponential convergence to some equilib- 
rium. However the rate of convergence tends to when a is close to ^ (in the 
case where J[fo] 7^ 0). So we do not expect to have a similar rate of convergence in 
the critical case. 

First of all, we know by Proposition l3.4l that the solution converges (in any H^iS)) 
to the uniform distribution as time goes to infinity. The goal of this section is to 
estimate the speed of convergence to this equilibrium. 

Proposition 4.3. Suppose that \\f{t) — IWh" is uniformly bounded on [to, +00) by 
a constant K, with s > '^^^2^^ ■ 

Then for all C > 1, there exists S > 0, such that if ||/(to) — 1||l2 ^ we have, 
for t ^ to, 

l|/(t)-l|U2^ ^ 



1 L 2(n^l) 

\/2(^ll/(to)-l|L2 "("+2) 

The constant 6 depends only on a, s, K , and C . 



it 



Proof. As in the previous section, we work on [to, +00). We write f = 1 + h and 
as in the previous case, we suppose that J[fo] 7^ 0. By the same argument used in 
Proposition 14.11 we have that J[f{t)] 7^ for all t > 0, so we define Q{t) as the unit 
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vector jTfjmjy- Similarly we denote (■) for the mean of a function on the unit sphere 
and cos for u ■ il. 

We have (h) = 0. Since Q is the direction of J[f] = ((1 + h)u) = (hu), we get 
that (hoj) = {hcos6)Q. 

We perform an expansion of the free energy and its dissipation in terms of h. We 
get, using (IHTj) and taking cr = ^, 

T{i + h) = - i{h') + j-^{h')) - lihcosor + omi)- 

Now we write a = n{h cos 9) and we define 

g = h - acos9 - |a^(cos^6' - i) - |a^(cos^6' - ^ cos 6'). (53) 

We have (cos"^ 9) = ^-^(^^^2) have used the formula (IH^ to compute = (cos^ 9)). 
Since we have (cos^ 6*) = {cos 9) = 0, and (cos^ 6') = ^, we get (g) = {gcos9) = 0. 
We will see that the terms of order 2 in g will not vanish in the expansion of the 
free energy and the dissipation term. But we will need to expand the free energy 
in a up to order 4, and the dissipation term up to order 6 in a. We have 

l{h') = + ^a' + ^^^a' + '^a'igcos' 9) + 0{a^g\U + (54) 
-W) = - 2n"(;+2) ^' - la^gcos" 9) + 0(||^||L + + a'lklU + c^'), 

We finally get 

-^(1 + h) = ^(/) + 4;^3(^a' + OdkllL + + «'ll^?lloo + a'). (55) 

Using the inequality a^h'^ ^ sa^ +(1 — s)f)~ for s G (0, 1), with a = a and h = ||5'||oo5 

2-|-— o I 1 

we get that a\\g\\l^ ^ |q;^ + lllfiflloo ' and a^lls-Hoo ^ + flls'lloo 

By Sobolev embedding and interpolation, as in the previous section, we have 

||^?||oo^C||(?||^, \\gyi , (56) 



with 1 - > f . 

Since a is controlled by using the definition ( l53l) of g, we have a bound 

2+- 

for ||5'||_H's on [tg, +oo), depending only on s and We finally get \\g\\oc ^ ^ C((7^)^, 
with > 1(2 + |)f > 1. 

So using fl54l) and (l55l) . we get that for any £ > 0, there exists 6 > such 
if \\h\\i2 ^ 6, we have 



(1 - e){{g') + ^ {h') ^ (1 + £)((/) + ia^) 
(1 - ^)i^{9') + 4^«') ^ -^(1 + h)^ ^^){2n'{n + 2)(/) + a'). (57) 

From that, up to take a smaller d, we obtain 

2nJ^(l + h) ^ (h^) ^ -i^ 2v/n(n + 2)J^(1 + h). (58) 
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We now estimate the dissipation term. We use the definition ( l32l) of !'(/) and the 
Poincare inequahty to get: 

V{f) = ((1 + h)\V{^Hl + h)- ((1 + h)u) ■ co)\') 
= ((1 + /i) I V(i ln(l + h)-{h cos e) cos 6*) |2) 

> ^(1 - l|/i||oo)(Nl + /i) - (ln(l + h)) - n{hcos9) cosO]^). (59) 
V ' 

S{h) 

We have 

= ln(l + h) - (ln(l + h)) - n{h cos cos 9 

= h-{h)-acose- - {h^)) + _ (/i3)) + ). 

We compute, 

h - (h) - acos6 = g + ^a^(cos^ ^ ~ «) + |"^(cos^ ^ - cos ^) 



2 

]{h' - (/i^)) = ia^ cos^' 9 + + c^Ml + «'ll^lloo + a^). 

So 

(^(/i)^) = ([^? + la%f^ - cos^)]^) + OiWgf + ahW + «1^lloo + a') 

= (9') + ;i3(^«' + 0{\\g\\l + a||^?||L + a%\\^ + «^). (60) 

As before, we get that ^ fa^ + l\\g\0 and a^g\\^ ^ fa^ + l\\g\0. 

Using (IMD, we get \\g\0 ^ C{gY, with > |(2 + |)f = 1. So using §9i) 
and ( I60l) . up to take a smaller 5, we have, for \\h\\L2 ^ 5, 

Vif) > (l - eV-^iig') + ,^a% 

Now for any C, C" > 0, if we take a and sufficiently small (so again up to take a 
smaller 6), we have that C{g'^) + ^ {C {g'^) + 0^)2. So we get 



.1 _ (/,2^) ^ _i ^ ^3 cos^)(cos2 ^ - i) + 0(||^f + a||(7|U + a' 



Putting this together with ( l57j) and the conservation relation (l42l) . we get that for 
any < e < 1, there exists 5o > such, as soon as ||/i||l2 ^ 5o, we have 

^^(/) = -!>(/) < - ^'"-/f--'' l-F(/)lt. 

(1 +£:)2^n(n + 2) 

Then we obtain, for all T such that \\h\\i2 ^ on [^Oi T], 

HfiT))-'^ ^ Hfm-^ + "l^^^^^M - to)- (61) 

(l+£)2 ^ji(n+2) 
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Then, using ( l58l) . we get that for t G [^O) 7"], 



n-e 



(l+e) 2y/ n{n+2) ' 



We write C = (a one-to-one correspondence between < £ < 1 and C > 1] 

and we get 



+ 



2(n-l) 



(t - to) 



So if we take 6 < min(<5n, — , 

on [to, T] for all T ^ to- Otherwise, the largest of such a T would satisfy 



(62) 



^o); and ||/i(to)||L2 ^ 6, we get that ||/i||l2 ^ 6q 



6o=\\h{T)\\L2^C 



< So- 



y/2{n+2)S 

So the inequality ( 162|) holds for all t G [to, +oo), which ends the proof. 



□ 



With this proposition, since / tends to the uniform distribution in any i/*(S), we 
get that for any r < , there exists to such that we have ||/(t) — 1||l2 ^ ^ ^ y 

for t ^ to. We can even get rid of the to in this inequality since for any r <r < ^(^^2) ' 
for t sufficiently large, we have , ^ 

As in the previous section, using interpolation to deal with the other Sobolev 
norms of the solution would lead, for any 77 > and t sufficiently large, to an 
inequality of the form ||/(t) — 1\\hp ^ Crit~2+'''i, But we can actually do slightly 
better. Indeed we have, following the notations of the proof and using (153|) . 



< |«||| COS 6*11 + 020^ + Cs\a\^ + Hfi'll^p- 

We have || cos9\\hp = {n — 1)2. We take to > satisfying the conditions of the 
proposition and such that ||/i||l2 ^ 5- We have that g is uniformly bounded in 
any if^(§), and so by interpolation, we have ||(?||_h's ^ for any 77 > 0. Now 

using (16T]) and (157|) . we get 



in}i{n+2) 



«')-5 ^ 



4(n-l)(l-£)? 



{l+e)^^n{n+2) 



which gives ||o||i2 = 0(t ^) and ^ — (i+e)^ n{n+2) — ^ finally, for any ?7 > 0, we 

2{n-l)(l-e)^(t-to) 



have that 



HP 



(l+£)^n(n+2) 



0(t ^^'^). This gives that there ex- 



(l+£)^n(n+2) 



2(n-l){l-e)3(t-to) 

ists ti ^ to such that for all t > ti, we have wp ^ (l+e)(n— 1)2 

' V 2(n-l){l-£)5(J-Jo) 

This is true for any e > 0. In conclusion, we have that for any r < n{n-i)p-'^ {n+2) ' 
there exists ti such that for t ^ ti, we have ||/(t) — 1\\hp ^ 
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4.4 Summary 

In summary we can state the following theorem: 

Theorem 4. Convergence to equilibrium. 

Suppose /o is a probability measure, belonging to (this is always the case 

for some s < — )■ 

Then there exists a unique weak solution f to Doi equation ([3]), satisfying the 
initial condition /(O) = /q. 

Furthermore, this is a classical solution, positive for all time t > 0, and belonging 
to C°°((0,+oo) X §). 

// J[/o] ^ 0, then we have the three following cases, depending on a. 

• If a > -, then f converges exponentially fast to the uniform distribution, with 
global rate {n — l){(r — ^) in any norm. 

More precisely, for all to > 0, there exists a constant C > depending only 
on to, s,p,n, and a, such that for all t ^ to, we have 

• If a < ^, then there exists G § such that f converges exponentially fast to 
Mf^Q^, with asymptotic rate roo{a) > in any norm. 

More precisely, for all r < roo{a), there exists to > (depending on fo) such 
that for all t > to, we have 

\\f{t)-M^n\\m^e-'^\ 

When a is close to ^ we have that To^io) ~ 2(n — 1)(^ — a). 

• If a = then f converges to the uniform distribution in any norm, with 
asymptotic rate (n+2) ^ 

More precisely, for all r < ^{n-iy-'^ {n+2) ' ^^^'^^ exists to > (depending on fo) 
such that for all t > to, we have 

\\f{t) -iWrn^ ^ 



rt 

If J[fo] = the equation reduces to the heat equation on the sphere, so f converges 
to the uniform distribution, exponentially with global rate 2na in any norm. 

For the subcritical case cr > ^, we used Theorem [21 In the case where p < — 
a simple embedding gives — l||_f/p ^ 11/(0 ^ 1|| -"-^ so we only have to show 

the result for p ^ We get 

The last inequality comes from the fact that /(to) is a probability density function, 
so /(to) — 1 is the orthogonal projection of /(to) on the space of mean- zero functions. 
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Using Proposition 12.61 we get 11/(^0)11/^^ ^ C'toll/o||if= in the case p ^ s. Otherwise 
we just use a simple embedding to get first ||/(to)||j?p ^ ||/(^o)||/i"'' and then by the 
same proposition ||/(to)||iyp ^ C'll/olli?''- 

Then the results in the case a < - and a = - are a summary of the conclusions 
of the two previous subsections. However, although it gives a clear understanding of 
how fast the solution converges to the equilibrium, in some sense, this summary is 
not as accurate as Propositions 14.21 and 14.31 which give a kind of stability: starting 
close to an equilibrium, the solution stays close. 

5 Conclusion 

In this paper, we have investigated all the possible dynamics in large time for the Doi- 
Onsager equation with dipolar potential. We have obtained a rate of convergence 
towards the equilibrium given any initial condition and any noise parameter cr > 0, 
for all dimension n ^ 2. 

The rate of convergence to the anisotropic steady state, in the case a < ^, 
depends on a Poincare constant which does not seem easy to estimate. A better 
knowledge of the behavior of this constant, for example as the noise parameter a 
tends to zero, would be useful to understand the limiting case a = 0, where we have 
existence and uniqueness of the solution. In this limit, the steady states are given by 
the sum of two antipodal Dirac masses (1 — a)6n + aS^Q with f2 G § and ^ a ^ |. 
We conjecture that if the initial condition is continuous (and with non zero initial 
momentum), then the solution converges to one of these steady states, with a = 0. 

It should also be possible to get the same kind of rates for the Maier-Saupe 
potential, but there the classification of the initial conditions leading to a given type 
of equilibria is much more difficult, in particular in the case where two types of 
equilibria are stable. 
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A Appendix 

A.l Using the spherical harmonics 

For the following we will use the spherical harmonics, so we recall some preliminaries 
results. We fix n ^ 2 and work on and its unit sphere S„_i. 

Definition A.l. A spherical harmonic of degree £ on Sn-i is the restriction to S„_i 
of a homogeneous polynomial of degree i in n variables (seen as a function — t- 
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which is an harmonic function (a function P such that AP = 0, where A is the 
usual Laplace operator in W^). We denote "H^"^ the set of spherical harmonics of 
degree £ on §n-i (including so they are vector spaces). 

We know that the space of homogeneous polynomials of degree i in n variables 
has dimension ("^1^^) (the number of n-tuples (ii,...z„,) of sum £). Writing an 

arbitrary homogeneous polynomial P of degree i under the form P = Yli=o Qi-i-^n^ 
with the polynomials Qi being homogeneous of degree i in the first n — 1 variables, 
and imposing that P is an harmonic function gives the following conditions (taking 
the term in X;-^), for i G |0, £ - 2l: AQe^i + {i + l){i + 2)Qf_i_2 = 0. Finally the 
polynomial P is only determined by the polynomials Qe and Qe-i in n — 1 variables, 
of respective degrees i and £ — 1. This gives the dimension of the space of spherical 
harmonics. 

Proposition A.l. The dimension ofH^/^^ is given by 

'^e - \ n~2 ) ^ \ n~2 ) - \ n~l ) \ n~l ) ' 

The second expression comes from two successive applications of Pascal's triangle 
rule, and will be useful in the following. It can also be seen by the following propertju: 
every homogeneous polynomial P of degree £ can be decomposed in a unique way 
as H+ \X\'^Q, where H is harmonic of degree £ and Q is homogeneous of degree £ — 2. 
Iterating this decomposition, we get 



P = Hi+ \X\^He^2 + \X\^H,_4 + 




even 
odd ' 



where the polynomials Hi are harmonic of degree i. This shows that any restriction of 
a polynomial on the sphere is equal to a sum of spherical harmonics (the terms 
are constant when restricted to the sphere). This gives, with the Stone- Weierstrass 
theorem, that the sum of spherical harmonics are dense in L^(S„_i) (since they are 
dense in the continuous functions). Together with the radial decomposition of the 
Laplacian A = -^dr{r"'~^dr) + ^A^ (where A^^ is the Laplace Beltrami operator 
on the sphere which is self-adjoint in L^(S„_i)), we get the following result: 

Proposition A. 2. The spaces 'Hf\ for£ e N, are the eigenspaces of the Laplace 
Beltrami operator A^ on the sphere for the eigenvalues —£{£ + n — 2). They 
are pairwise orthogonal and complete in L^(§„„i). 

We can construct a basis of 'H}f^ by induction on the dimension, using the sepa- 
ration of variables. We describe this construction and will use it in the following. 



^This can be shown using the appropriate inner product (P, Q) i— >■ P(D)Q on the space of 

homogeneous polynomials P of degree where P(D) is defined as „„„ if P = V . -X"", 

■■■"x„ 

and extended by linearity (so for example, we have that |Xp(D) — A). If we denote by E the space 
of polynomials of the form P = jXpQ, with Q of degree i ~ 2, then the orthogonal of E consists 
in all the polynomials P such that for ah Q of degree £-2, we have (|XpQ)(D)P = Q(D)AP = 0, 
that is to say in all the polynomials P such that AP = 0. So the claimed decomposition is just 
the orthogonal decomposition, on E and E^. 
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For a given unit vector e„ G M", we take an orthonormal basis (ei, . . . , e„) of M". 
Any u G Sra-i \ {Gni—e-n} can be written a; = cos^e„ + sin^t), with ^ G (0,7r) 
and V G S„-2- We identify M"~^ with the vector space spanned by (ei, . . . ,e,i_i). 
The special case n = 2 works if we consider 5*0 = {ci, — Ci}. 

By convention, the only spherical harmonics on Sq are the constant functions (of 
degree 0) and the functions ci H- c, — ei i— ?• — c (of degree 1). 

Now, for n ^ 1, we choose an orthonormal basis (^^, . . . ^) of 7im~^^ for 

any m G N and we have the following result: 

Proposition A. 3. There exists polynomials Qe,m of degree i — m such that if we de- 
note y,^(u;) = g,,^(cos^) sin'^eZ^iv), then the Y^^ for m G |0,£], G |1, fc^^"'^] 
form an orthonormal basis of H^^^ . 

Proof. Writing Y^^{u) = Qi^rn{cos9) sin"^ 9Z!^{v) and asking it to be a spherical 
harmonic is equivalent to the following linear ODE for Qf.m (we recall that the 
Laplace-Beltrami operator is given by sm'^~"' 9d0{sin^~'^ Ode) + i^^e^i^ this coor- 
dinates): 

sin^-" dei - sin"+'"-^ OQ'g^^{cos 6) + m cos 9 sin"+'"-3 6'Q£,„(cos 9)) 

- m(m + n- 3)Q^,„(cos 9) sin'"-^ 9 = + n- 2)g^,^(cos 9) sin™ 9. 

We write x = cos 9 and this equation transforms into 

(1 - x'^)Ql^ - (n + 2m - l)xQ'^„^ + {i - m){i + n + m - 2)Q^,„ = 0. 

This equation is a particular form of the Jacobi differential equation, where the two 
parameters a and /3 are equal (also called Gegenbauer differential equation). One 
solution of this differential equation is a polynomial, called ultraspherical polynomial 
(a particular case of the Jacobi Polynomials, also called Gegenbauer polynomials), 
and denoted Pj^^^ following the notation of Szego in [T7]. Precisely, it satisfies the 
differential equation 

(1 - x'^)y" - (2A + l)xy' + t{i + 2\)y = 0. 

Taking A = m — 1 + | and i = i — m, we get a solution Qe^m = C(e,mPi-m ^ > 
where a^^^ is a positive constant of normalization, such that is of norm 1 

in L^(S„_i). We have to be careful here because Pj^^'' is not defined for A = 0, and 
so the only special case is n = 2, m = 0, for which we have a solution Q^ q = \/2T£, 
where r£(cos^) = cos £9 (the Chebyshev polynomial of first order of degree i). 

So for a fixed £, we have constructed a family of spherical harmonics Yj^^ of 

degree i for m G |0,£],A; G |l,A;m ^^J- They are pairwise orthogonal in L^(S„_i) 
since the are pairwise orthogonal in I/^(S„_2). The size of this family is exactly 

m=0 m=0 

which is the dimension of Hf^\ so we get that the Y^^ for m G |0, A; G |1, A;^^ 
form an orthonormal basis of H^^\ □ 
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From now on, we will use the construction done in the proof. We have that, for a 
fixed m ^ 0, the polynomials Qe^m for i ^ m are a family of orthogonal polynomials 

1 . n— 1 

for the inner product (P, Q) ^ P{x)Q{x){l - x^r-'^~dx. 

We will use three properties on the Gegenbauer polynomials (see [TT]) for the 
following, for z ^ 0, A 7^ 0, and A > — | (with the convention =0): 



^(p/^)(.))^(i - xy-.dx ^ (64) 

{^ + = 2(. + A)XPW - (. + 2A - 1)pS (65) 

(1 - X'){Pi'^y = ^^-1^ ((z + 2A - l)(z + 2A)p£i - z(z + 1)P,51) (66) 

We have the following normalization for the Qe^m'- 

' Ql^{x){l - x^)— i+Vdx = f\l- x^) V-Mx. 
1 ' J -I 

This gives the following relation, together with ( !64|l : 

2 _ (£+t)(f+l-m) 2 

By the previous construction, we can decompose g = Ylk i m'^XrrXt,m '^^ 
have Jg ~ Ylikim k^,mP- Since g is of mean zero, we have cjg = (the only 
spherical harmonic of degree is the constant function 1). So from now, the in- 
dices fc, m of the sum will mean £ > 0, m G |0, £], A; G |1, km 

We decompose in the same way h = Ylkim'^imXi'm- S^^^ ^ fii'st formula, in 
the form of a lemma. 

Lemma 3. We have 

1 



where = ^^^=^1^+7^" ^ ^ 1. 



Proof. We have 

en ■ VY^,m = - sin = [(1 - X^)Q',^^ - mXQ,,^] (cos 9) sin™ ^Z^(z;), 

and using the inductions formulas (lUUj) and f lU7|) . we get 

(1 - X^)(5^„j - mXQi^rn = ^[be-l,ra{i + n - 2)(5^_i - bi^rniQe+l,m], (69) 

where 6^,^ is given in the statement of the lemma. In the special case n = 2, m = 0, 
using the formula Q£^q{cos6) = cos £6 gives the same formula as (169|) . with = 1- 
So we have that Jg ^ e„ ■ VY^^^^Y^^^/ can be non-zero only ii m = m' , k = k' , 
and £ = ^' ± 1. By bilinearity, together with the fact that form an orthonormal 
basis, this gives the claimed formula. □ 
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Now we have all the tools to prove Lemma [U (we recall it here) . 

Lemma 1. Estimates on the sphere. 

1. If h in H~^~^^{E>) and g in H^(S>), the following integral is well defined and we 
have 



gVh 



H- 



where the constant depends only on s and n. 
2. We have the following estimation, for any g G H^^^{Ei) 



[ gVi-AYg 



^C\\g\ 



where the constant depends only on s and n. 
3. We have the following identity, for any g G H~ 

gVA-\g = 



n-3 
2 , 



Proof. Using Lemma [3l we get 



^- J2 

k,£,m 
1 

+ 2 



i+1 \ A 



E 

k,£,m 



i+n-2 \ X. 



-s + l 



^C\\9\ 



h\ 



where Xi = i{£ + n — 2) (the eigenvalue of —A for the spherical harmonics of 
degree £). The last line comes from the fact that the sequences ■ / ^+^~'^ ( h+i. 



and 



e+n-2 \ Xi+i I bounded (they tend to 1), together with a Cauchy-Schwarz 
inequality. This gives the first part of the lemma, since this is true for any unit 
vector e„. 

Now we take h = {—A)^g, which is replacing d'l^ by A^c^^ in Lemma [3l We get 
en ■ [ g^i-AYg = Yl \hn.4+i,^clJ{l + n- l)A,%i - £A|] 



2 „fc 
t '-'0 



\lclJ\{£ + n-l) 



Indeed we have that = 1 — | 



+ 0(^2), SO |(£ + n- 1) 



IS 



Af t V/' IV ' y \i J \A^_|_i 

bounded (it tends to (n — 1) + 2s). Since this computation is now valid for any unit 
vector e„, this gives the second part of the lemma. 
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The last part is straightforward by taking h = ^n-id with Lemma [3l According 

k ^ 1 

i,m £(^e+l)...(e+n-2) 



to the definition given in ([8]), we have d'^^= ]f,„_2) 4m- We get 



"^^"-1 k,l,m 

which is true for any unit vector e„. □ 
A. 2 Analyticity of the solution 

Following |3], we will show that the solution belongs to a special Gevrey class. We de- 

fine the space Gr, as the set of functions g (with mean zero) such that A^^-^^e^^~^^^ g 
IS m L2(§). Using the notat ions of the previous proof, this is an Hilbert space 
associated to the inner product 



The norm on this Hilbert space will be written || ■ Wg^- 

We take r a function of t, we will denote its time derivative by f. For a given 

solution / = 1 + ^, we put /i = A~\e^'''^-^'^^ g in ffT3l) . 
The left-hand side is 



E 2ry/i{e+n~2) i 
£(£+l)...(^ + n-2)^^''"dt'^^''" 



Id g2rv'<?{f+n-2) f k \2\ ■ g2rv^£(f +n-2) | ^ |2 

9d7^^ 



Using Lemma [HI we get 



11 2r^(£+l)(£+n-l) _ 2ry/e{l+n~2) 

e ■ (o VA'^ p2r(-A)^ \ _ tsr^ 1^ k k ^ ^ 

en {g,\/^n-ie 9) - ^ 2^ b,^mC,+, r^C,^^ (£+ 1) + ^ _ 2) 

k,£,m 



k,l,m 

X 



4/ £(^+n-l)~ / r(^(^+l)(^+™-l)-^^(£+n-2)) _ ^-r(^(^+l){^+n-l)-^£(^+n-2)) 

V (t+m+n-2) ' t ^ 

^sinh(r(\/2n- Vn- 1))||(-A)^^||g^. 



Indeed the expression ^/(^ + 1)(£ + n — 1) — + n — 2) is a decreasing function 
of £ ^ 0. Since this is valid for any unit vector e^, we get 



JM-((?,VA;\e2^(-^)^(7) ^sinh(r(v^-v/;^))||(-A)i^? 



i2 
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Now since ||(-A)i(?|||^ ^ _i=||(_A)l^||2,^, and | J[/i]| ^ J [^]|, we finally 

get 

Ij^Mk + k - + sinh(r(v^ - V^^)))] ||(-A)5^|||^ ^ 

As soon as r + sinh(r(-\/2n — y/n — 1)) < [a — e)^/n — 1 and r is bounded in time, 
we have that Hq'IIg^ is uniformly bounded (and we can indeed take r{t) = 5min(l, t) 
for 6 sufficiently small), provided go is in Gr{o)- So if we have r(0) = 0, we only need go 
to be in H ~ . If it is not the case, by instantaneous regularization (Proposition [22]) 
we have it for any time t > 0. Since Gr, for r > 0, is a subset of the set of analytical 
functions on the sphere, we get that any solution becomes instantaneously analytic 
in space. 
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